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Abstract 

The photon sphtting 7 ^ 77 in a time-independent and inhomogeneous mag- 
netized background is considered when neutral and ultra-hght spin-0 particles are 
coupled to two-photons. Depending on the inhomogeneity scale of the external field, 
resonant photon splitting can occur. If an optical laser crosses a magnetic field of 
few Tesla with typical inhomogeneity scale of the order of the meter, a potentially 
observable rate of photon splittings is expected for the PVLAS range of couplings 
and masses. 
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Ultra-Light scalar /pseudoscalar particles have escaped, so far, direct detection. There 
are, however, theoretical reasons to believe in their possible existence. For instance, 
they are naturally predicted by diverse theories beyond the standard model endowed 
with spontaneously broken (global) continuous symmetries. Because of the lack of any 
direct detection, these particles must be very weakly coupled to standard matter fields. 
A pivotal example along this direction is the (pseudo-scalar) Nambu-Goldstone boson 
associated with the spontaneously broken Peccei-Quinn symmetry [1]. The mass of this 
particle, customarily called axion [2], is believed to be in the meV range [3,5]. 

As far as the electromagnetic interactions are concerned, the Lagrangian densities for 
an ultra- light pseudo-scalar {ipp) or scalar {^ps) field can be parametrized, for the purposes 
of this discussion, by two different couplings, i.e. A and A: 

Lp^-^F^J^^^P, Ls^-^F,,F^^^S: (1) 

where F^^, and F^^ are, respectively, the Maxwell field strength and its dual. 

Astrophysical constraints [5] demand that the axionic coupling A should be of the 
order of 10^° — 10^^ GeV, implying, together with the smallness of the corresponding 
mass, that the axion is (almost) stable in comparison with the age of the Universe. The 
axion would then be a potential dark matter candidate. 

Recently, the PVLAS experiment [6], has reported the first evidence of a rotation of 
the polarization plane of light propagating through a magnetic field. According to the 
standard treatment [7], these results would imply, if confirmed, the presence of a very light 
pseudoscalar particle (axion) whose inferred mass m and coupling A will be, respectively, 
C(meV) and 0(10^ GeV). The purpose of the present paper is to test this claim by 
considering a complementary effect that has not received, so far, specific attention. Our 
logic can be summarized, in short, by Fig. 1. Suppose that a laser beam passes through 
an inhomogeneous magnetic field. Now, if the magnetic field can absorb momentum in 
a continuous manner, the (off-shell) axions can be reconverted into photons. In fact, 
the two-photon couplings in Eq. (1), could also induce another indirect effect which 
is the photon splitting 7 — 77 in an external (time-independent) magnetic field. In 
a QED framework [8], the rate of photon splitting is suppressed by (B/B^™)®, where 
B§^^ = ml/e ~ 4.4 x 10^ Tesla. Owing to the largeness of B^™ the resulting effect is 
extremely minute for typical laboratory (i.e. O(Tesla)) magnetic fields. The question we 
ought to address is therefore rather simple: how many photon splitting events are expected 
in the situation described by Fig. 1? If the magnetic field would be completely static 
and homogeneous the answer to this question would be only academic since, due to the 
translational invariance of the full system, the process could only proceed by taking into 
account photon dispersion effects [9]. Besides that, the process would also be suppressed 
by two extra powers of A. 

If, on the contrary, the (classical) magnetic field is inhomogeneous, it is plausible to 
expect that the momentum of the photon beam could be partially absorbed in a contin- 
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Figure 1: The Feynman diagram describing the contribution to the photon sphtting amphtude 
in an external magnetic field. The spin-0 particle {(j)) is coupled to two-photons as in Eq. (1). 
The dashed horizontal line indicates the propagator of the spin-0 particle, and the vertical dot- 
dashed line stands for the external magnetic field insertion. The bubble represents the magnetic 
form factor; k, q, and ki^2 denote, respectively, the initial, intermediate and final four-momenta. 

uous manner. In the latter case the answer to the aforementioned question is different 
since possible resonant effects must be taken into account. Indeed there are two relevant 
physical scales in the problem at hand: L (i.e. the magnetic inhomogeneity scale of the ex- 
ternal background field) and p (i.e. the momentum absorbed by the magnetic field at the 
resonance). If L ~ 0{l/p) the photon splitting production rate is said to be resonantly 
amplified. Amusingly enough the typical length scale L turns out to be of the order of the 
meter for the PVLAS range of masses [6] and taking, as incident beam, an optical laser. 
While inhomogeneous magnetic fields have been considered in order to produce real (i.e. 
on-shell) axions by Primakoff effect [3, 10] (see also [11, 12]), their possible relevance for 
photon splitting has not been taken into account so far. One of the results of the present 
paper is that the PVLAS axions should also produce an observable excess of photon- 
splitting events in comparison with the (minute) QED background previously mentioned. 
Consider therefore the process illustrated in Fig. 1, i.e. jlk) — > 7(^1) 7(^2) ? where k and 
ki^2 are, respectively, the four momenta of the initial and of the final photons. Suppose 
also that the photons exchange momentum with the external magnetic field (represented 
by a cross in Fig.l) which is assumed to be static (i.e. time- independent) but spatially 
inhomogeneous. Consequently, while the total energy of the reaction is conserved (owing 
to the time- independence of the external field), the three-momentum can be absorbed by 
the (inhomogeneous) external field. 

Defining the four momenta of the final photons as ki = {ui, ki) and k2 = {002, k2), the 
differential cross section for the scattering of a photon by a time-independent external 
field is: 



where M(7 77) is the matrix element of the process and the 1/2-factor arises since 



dcr= - H |M(7 ^ 77)1' (2vr)5(^ - - ^2) 
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the photons of the final state are indistinguishable. According to Fig. 1 the scattering 
amplitude for 7 — >■ 77 in a magnetic field can be written as 

sm ^ % 
Mij^jj) = —.B(p,j^-^^^^^-^^^,^V^(M,k.), (3) 

where m and F are, respectively, the mass and total width of the exchanged boson; 
moreover 6 is the angle between the direction of the external magnetic field B and the 
incoming photon three-momentum k. The term B{p) is the Fourier transform of the 
projection of B(a:) along the polarization vector of the incoming photon, namely B(p) = 
/ \n ■ B(a;)| e^^'^cPx, where p = k — ki — k2 is the momentum absorbed by the external 
field and ft is the unit vector parallel to the direction of the incoming polarization. In 
Eq.(3), the term 

V^^ik,, k2) = J el k^ el) e^,^p , (4) 

represents the boson-77 vertex contribution where Si{ki) and £2(^2) are the polarization 
vectors of the two final photons and e^ij^p is the total antisymmetric tensor. 
Inserting Eqs.(3) and (4) into Eq.(2) the differential cross-section becomes 

8A4((A;i + A;2)2-m2)2 + F2m2' 

X (2'K)d(u; — LOi — 102) ; ; r^. (5) 

As previously mentioned, for the reaction to proceed, the three-momentum p must be 
absorbed by the external magnetic field. In fact, if the external field is fully homogeneous 
not only in time but also in space, we will have that, in Fourier space, |B(p)| oc 6^^\k — 
ki — k2). This occurrence demands that Eq.(5) is proportional to {ki + k2y 5^'^^k^kl — k2)■ 
Then the total cross section vanishes after integrating over the total phase space. 

Naively, if the magnetic field has a finite extension of order L and vanishes outside 
of it, the translational invariance is broken and the external field could easily absorb 
momentum of magnitude ~ 1/L. To model this situation, consider, as an example, a 
magnetic field configuration with Cartesian components B{x) — (B-r(x), 0, Bz{x)), where 



B,(f) = B exp f-^ - , B.(f) = ^ £^ dxB,{x) . (6) 



It is easy to check that B{x) is indeed solenoidal. By choosing the three-momentum k 
of the incoming photon along the x-axis, only the z component of the magnetic field will 
enter Eq. (5) where |B(p)| \J d^xBz{x)\. 

Let us then compute our main observable, i.e. the number of photon splitting events 
taking place inside the magnetic field. This quantity, denoted in what follows by N^^, 
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is simply the product of the cross section o (obtainable from Eq.(5)) times the flux of 
incoming photons $ = N^{t)/A , i.e. N^^{t) = a^. The quantity N^{t) is the number of 
incoming photons per unit time crossing the magnetic field and A — nB? is the surface 
spanned by the magnetic field on the the plane orthogonal to the direction of the incoming 
photon momentum (in our frame the y — z plane). Inserting now Eq.(6) into Eq.(5), 
the differential number of photon splitting events per unit time is given by dNjj{t) — 
N^{t) ciP(7 77), where 

c^P(7-77) = ((fc, + kly-m^y + Pm^ " I 1 >< 




(7) 



2a;i(27r)3 2u;2(27r)3 ' 

is the differential probability of 7 going to 77. In Eq. (7), the absorbed momentum is 
p = u—kix—k2x- The quantity Ai(p) is nothing but Ai(p) = {L/y/n) exp (—L^ p^), which 
is a well known representation of the Dirac-delta function 6{p) in the limit lim^^oo Al{p) = 
S{p). To simplify the problem, the shape of the magnetic field can be chosen in such a 
way that R 00 while L stays finite: in this way the total momentum will be absorbed 
only along the initial beam direction x. In practice this means that L, i.e. R much 
larger than the characteristic inhomogeneity scale along the x direction. 

After using the Dirac-delta functions for the partial integrations a change of variables 
allows to write the total probability as 

^ = 25^^^^^^'^^+^'^'^^^ ' 
where the following integrals have been introduced 

x,H=./^d.,/^ (10) 

F(m') = L^e-''''''\ p^u,(l-^l-^-'ia'). (11) 

The relevant physical regime for the present purposes is the one where uj > m. This limit 
is realized, for instance, when the incident photon beam is in the optical and the scalar 
mass in the PVLAS range. If m/uj <^ 1, the Breit-Wigner distribution can be easily 
integrated in the thin-resonance limit (since F = m^/(647rA^) -C m), and, from Eq. (8) 
the result is simply 

P(7 - = ^ ^(^) BR(V'P - 70) (1 + 0(I) ) . (12) 
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(13) 



where BR((/9p — ^ 77) = m^/(647rA^r) is the branching ratio of the pseudo-scalar decay 
in 77 and terms of order 0{m'^/uj'^) have been neglected. Furthermore, in Eq. (13) p is 
the absorbed momentum at the resonance, i.e. 



Note that J-'{x), being dimensionless and normalized to 1, measures the suppression of 
P{l ^ 11) when the magnetic inhomogeneity scale L is larger than 1/p. The limit 
a; ^ m is naturally implemented for the case of PVLAS axions when uj is the energy of 
the (optical) photon beam, and in this case p ~ w? / {2uj). Notice that the leading term of 
the integrals in Eq.(8) is proportional to l/F, enhancing the suppression induced by the 
1/A^ coupling. In particular, when the pseudo-scalar has only the decay channel in two- 
photon, BR((/7p — > 77) = 1, Eq. (12) reproduces the probability of photon conversion 
in on-shell axions [3]. Such an occurrence is a direct consequence of the Breit-Wigner 
distribution in Eq.(5). It is relevant to mention that if a; < m (which is opposite to the 
limit where Eq. (12) has been derived), the region of the pole in the integrand of Eqs. 
(9) and (10) does not contribute to the total integral. Consequently, the total probability 
is simply suppressed as B^/A"^. 

Generalizations of the result (12) to the scalar case are possible. In the scalar case the 
boson-77 vertex of Eq. (4) is different and the coupling will be dictated by A. Conse- 
quently, the relevant geometric set-up will be the one where the incoming photon polar- 
ization is orthogonal (and not parallel) to the orientation of the external magnetic field. 
Moreover, the dependence of the branching ratio for the decay in 77 upon the mass, 
couplings and width of the scalar particle has the same analytical expression as in the 
pseudo-scalar case. With these caveats the numerical results are the same if A = A. 

According to Eq. (13), the resonant effect in Eqs. (12) is achieved when x = 
corresponding to the maximum of J^{x) as a function of x- the case of an incident 
optical laser with a; ~ 1 eV, and for a magnetic field of 0{1 Tesla), Eq. (12) implies that 
P{l — 77) ~ 4 X 10~^^ for mass and coupling in the PVLAS range [6], i.e. m ~ 10~^ 
eV and A = 10^ GeV. For these figures, the characteristic magnetic length L will be 
macroscopic, i.e. L ~ 0.5 m. For a Nd:Yag laser with average power of 1 Watt (and typical 
wavelength of A ~ //m), we will have that N^{t) ~ 10^^/s. Consequently, the number of 
photon splittings per second will be N^^ ^ 4 x 10^ which is potentially observable. As it 
can be appreciated from Fig. 2 the PVLAS region can be excluded, at 95 % C.L. (10^ sec 
of integration time), for magnetic inhomogeneity scales smaller than about a meter. 

Different assumptions on the magnetic field profiles lead, up to numerical factors 
(9(1), to the same results. This has been cross-checked with a variety of profiles where the 
magnetic field vanishes at infinity. Moreover, the different shapes of the magnetic field will 
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Figure 2: Exclusion regions for photon splitting (above the lines) at 95% C.L. in the (m , A~^) 
plane, for an optical beam of frequency uj = 1 eV, and a magnetic field B = 1 T of different 
length (L), namely L/m=0.5, 1, 2, and 5. The point with error bars corresponds to the PVLAS 
range. The integrated time is 10^ sec. 



modify the form of J-'{x) and will therefore change the suppression of the probability when 
the magnetic inhomogeneity scale is larger than the inverse of the resonant momentum p. 

It is instructive to discuss briefly the massless limit of Eq.(12). In this limit, as a 
consequence of the Breit-Wigner distribution (see e.g. Eq.(8) ) the resonant pole is absent 
and the photon splitting probability is discontinuous in the limit m —>■ 0. In the massless 
case, the total probability of photon splitting should then be calculated by setting m = 
inside the integrand of Eq. (8) and the final result can be expressed as 

- ■ (15) 





256 7rA4 

4a;(l— x) 

where ^ = uj L. If, as previously assumed, u; = 1 eV and L ~ 0(m), we have ^ 1, and, 
in this limit, T2,{C) — ^ V^27r. This occurrence implies then 



which is very suppressed due to the fact that the resonant effect has disappeared for m = 0. 
As a useful comparison we can also report the probability of conversion of a photon into a 
(massless) pseudo-scalar which turns out to be P(7 ipp) = B^L^ 7r/(8A^). This result 
follows from the same external magnetized background assumed in Eq. (6). 

Of course, the aforementioned discontinuity is present if dispersion effects, associated 
with the self-energy corrections due to the interaction with the external field, are ne- 
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glected. These effects typically cause a shift of the squared mass, i.e. = 
(m^ + \Jm'^ + 4m|)/2, where uib = ^JBcj/A. Corrections, of order O(mB), arise as well 
for the photon self-energy. If the mass is in the PVLAS range (mpvLAS ~ 10~^eV) for 
Tesla magnetic fields and couphngs A ~ 10^ GeV, uib — 10~^ eV <^ mpvLAS a-nd disper- 
sion effects are negligible for resonant photon splitting. 

As already mentioned PVLAS results lack independent confirmations. In the present 
paper it has been shown that, for the PVLAS range of masses and couplings, an observable 
rate of photon splittings is expected when the external magnetic field can absorb three- 
momentum in a continuous way. Therefore, to confirm, in an independent channel, the 
PVLAS findings, resonant photon splitting represents an intriguing option. 
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